Contribution of the massive photon decay channel to neutrino cooling 

of neutron stars 



D.N. Voskresensky] 

Moscow Institute for Physics and Engineering, Kashirskoe shosse 31, 115409 Moscow, Russia 
GSI, Planckstr. 1, D-64291, Darmstadt, Germany 

E.E. Kolomeitse'v[j], B. Kampfer 
Institut fur Kern und Hadronenphysik, Forschungszentrum Rossendorf e. V., 
PF 510119, D-01314 Dresden, Germany 

Institut fur Theoretische Physik, TU Dresden, 
Mommsenstr. 13, D-01062 Dresden, Germany 



We consider massive photon decay reactions via intermediate states of electron-electron-holes and 
proton-proton-holes into neutrino-antineutrino pairs in the course of neutron star cooling. These 
reactions may become operative in hot neutron stars in the region of proton pairing where the 
photon due to the Higgs-Meissner effect acquires an effective mass m 7 that is small compared to 
the corresponding plasma frequency. The contribution of these reactions to neutrino emissivity is 
calculated; it varies with the temperature and the photon mass as T^m^e-™-^ for T < 
Estimates show that these processes appear as extra efficient cooling channels of neutron stars at 
temperatures T ~ (10 9 - 10 10 ) K. 

PACS number(s): 21.65.+f, 95.85.R, 97.60.1 

key words: nuclear matter, neutron stars, neutrino cooling 



* e-mail: voskre@rzri6f.gsi.de 
t e-mail: kolomei@tpri6f.gsi.de 



1 



I. INTRODUCTION 



The EINSTEIN, EXOSAT and ROSAT observatories measured surface temperatures of certain neutron stars and put 
upper limits on the surface temperatures of others ( cf. Ref . and further references therein) . Data on the supernova 
remnants in 3C58, the Crab, and RCW103 indicate rather slow cooling, while the data for Vela, PSR 2334+61, 
PSR 0656+14, and Geminga point to significantly more rapid cooling. In the so-called standard scenario of neutron 
star cooling, the most important channel up to temperatures T < (10 8 — 10 9 )K corresponds to the modified URCA 
process nn — > npev. Rough estimates of its emissivity were first made in Ref. (2| . Friman and Maxwell in Ref. (^] 
recalculated the emissivity of this process in a model, in which the nucleon-nucleon interaction is treated with the 
help of slightly modified free one-pion exchange. Their result for the emissivity, £„ M , proved to be an order of 
magnitude higher than previously obtained. The value £„ M was used in various computer simulations resulting in 
the standard cooling scenario; see Ref. Q for example. Subsequent work (^-Q] took in-medium effects into account in 
A^TV-interaction, showing that the emissivity of the modified URCA process depends heavily on neutron star mass. 
For stars of more than one solar mass, the resulting emissivities turned out to be substantially higher than the values 
given by £„ M ■ 

These and other in-medium effects were recently incorporated in the computer code in Ref. || leading to a new 
scenario of neutron star cooling. For low-mass stars numerical results of the new and standard scenarios more or less 
coincide. In the present work, we continue to look for enhanced reaction channels. To demonstrate the efficiency 
of new reaction channels, we compare the results with the emissivity £„ M , which dominates cooling in the standard 
scenario over temperature range under consideration. 

Besides the modified URCA process, the standard scenario numerical codes also include neutron and proton 
bremsstrahlung processes nn — > nnvD and np — > npvv, which in all models lead to a somewhat smaller contri- 
bution to the emissivity than the modified URCA process §,|J|§. Also included are processes that contribute to 
emissivity in the neutron star crust. These are plasmon decay 7 p / — > v v |io|Jll| , electron bremsstrahlung on nu- 
clei -> e Avv [ |12 | Jl3|Jl~l| ] , electron-positron annihilation e e + — > vv pl|J15[ |, and photon absorption by electrons 
7e — » e vv ]l6| , |ll|[i7| |~Numerical simulations show that the latter two processes contribute only negligibly to the crust 
neutrino emissivity at the temperatures under discussion in this paper and they always contribute negligibly to the 
full neutron star's emissivity; see Fig. 7 of Ref. pl]| . 

When the temperature decreases, it is energetically favorable for neutrons to pair in the neutron star interior and 
inner crust and for the protons to pair in the star's interior. In a system with nucleon pairing the emissivity of 
the modified URCA process is suppressed by a factor exp(— (A„ + A p )/T) ||, where A„ and A p are the respective 
neutron and proton gaps, defined by Aj(T) = Aj(0) (T c ^ — T) T~ t 9(T C) j— T) (here Q(x) is the Heaviside step function, 
i = {p, n}, and T c ,i is the corresponding critical temperature for nucleon pairing). At temperatures T -C T CjP ,T c , n 
the process becomes marginal. Nevertheless, this star's interior process still dominates those of crust cooling up to 
temperatures T ~ (10 8 - 10 9 ) K, depending on the values of the gaps; see Fig. 7 of Ref. O. For T < (1 - 3) • 10 s K 
cooling in the standard scenario is largely dominated by the photon emission from the neutron star surface. 

In the present work we look for more efficient cooling processes at T < T cp ,T cn . We analyze photon decay 
into neutrino-antineutrino pairs. The related processes "/e — ► e vv and 7p — > pvv turn out to be suppressed by 
several orders of magnitude compared to those under discussion, due to the lack of free final states in degenerate 
fermionic systems, and are therefore not considered here. The contribution of photon decay via electron-electron-hole 
intermediate states for the case of a normal electron plasma in white dwarfs and neutron star crusts has been calculated 
by several authors (see Ref. Jl0[ for further references). In an ultrarelativistic electron plasma, a photon acquires an 
effective in-medium plasmon dispersion law with a gap equal to the electron plasma frequency uj p i ~ 2e p e /^/3n, 
where e is the electron charge and p e denotes the electron chemical potential (we employ units with K = c = 1). 
Therefore, the contribution to the emissivity of the cited process is suppressed by a factor exp(— u> p i/T). Nevertheless, 
in white dwarfs and neutron star crusts, the electron density is not too high, and the process is still effective. In 
neutron star interiors, the electron density p e is equal to the proton density p p by virtue electrical neutrality, and 
along with (3 stability one obtains a relation for the total density 

Pe = p p ~ 0.016 po (j^j , (1) 

where po — 0.17 fm~ 3 denotes the nuclear saturation density, and we use the values of the neutron and proton Fermi 
momenta S, pp n — 340(p/po) 1 ^ 3 MeV and pp p = // e ~ 85(p/ po) 2 ^ 3 MeV. Thus, at typical densities for neutron 
star interiors p ~ po, the value of the electron plasma frequency is high, e.g., oj p i(pq) w 4.7 MeV for p ~ po, and 
at temperatures T < T Ctn ,T CtP < ui p i the process 7 Wp , — > ee _1 — > vv, where the superscript —1 denotes the hole, is 
strongly suppressed. We therefore seek another process that can contribute to rapid cooling. 
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We exploit the fact that, contrary to a normal electron plasma, in superconducting proton matter, due to the 
Higgs-Meissner effect, the photon acquires an effective mass that is small compared to the plasmon frequency. In 
the region of proton pairing at T < T c , p , we therefore find that new decay processes of massive photons (-f m ) 
via electron-electron-hole (ee _1 ) and proton-proton- hole (pp~ l ) intermediate states to neutrino-antineutrino pairs, 
7m — * ee _1 +pp~ x — > WW, I — {e,/i,r}, can dominate neutron star cooling at certain temperatures. These processes 
are determined by the diagrams 




In the first diagram, the solid lines in the loop are related to Green's functions of non-superfluid relativistic electrons. 
In the second and third diagrams, the solid lines in the loops correspond to superconducting nonrelativistic protons. 
The distinct orientations of arrows indicate that the second diagram is calculated with so-called "normal" Green's 

functions — * , which become the usual Green's functions for normal Fermi liquids in the limit A p — * 0. In contrast, 

the third diagram is built up with the "anomalous" Green's functions ■< >-and — »— * — , which are proportional to 

the proton gap. Therefore the contribution of the third diagram vanishes for A p — > 0. The fat vertices in the second 
and third nucleon diagrams include nucleon-nucleon correlations. 

The contribution to neutrino production matrix elements of the third diagram and terms proportional to the gap in 
the second diagram is as small as (A p /e^ p ) 2 <C 1 for T < T CyV <C cf p (here tFp is the proton Fermi energy), compared 
to the contribution of the second diagram calculated with Green's functions of the normal Fermi liquid. To this same 
accuracy, we drop the third diagram and use the Green's functions of protons for the normal Fermi liquid^] in the 
second diagram. We thus calculate the emissivity according to the first two diagrams, assuming A p = in the second 
diagram but taking into account that the photon dispersion relation is changed due to proton superconductivity. 

Our paper is organized as follows. In Sec. 2 we show that in the region of proton superconductivity due to the 
Higgs-Meissner effect, the photon spectrum is rearranged, and instead of the plasmon gap the photon acquires a mass, 
which is now determined by the density of paired protons. In Sees. 3 and 4 we demonstrate the efficiency of these new 
processes in the course of neutron star cooling. The emissivity corresponding to the above diagrams is calculated and 
compared with the emissivity of the standard URCA process and photon emissivity from the neutron star surface. In 
Sec. 5 we detail our conclusions. 



II. PHOTON SPECTRUM IN THE SUPERCONDUCTING PHASE 

As is well known fll8|| , the photon spectrum in superconducting matter and in a normal plasma are substantially 
different. In the superconducting matter considered here we deal with two subsystems. The normal subsystem 
contains electrons and non-paired protons and neutrons, which are present to some extend at finite temperatures. The 
superfluid subsystem contains paired protons and neutrons. In the presence of a superconducting proton phase, normal 
currents associated with both electrons and residual non-paired protons are fully compensated by the corresponding 
response of the superconducting current Jl8|,^0|,|l| ; otherwise there would be no superconductivity. What remains after 
this compensation is a part of the superconducting current. The resulting photon spectrum is thereby determined by 
the inverse of the London penetration depth ( due to the Higgs-Meissner effect Ul), but not by the plasma frequency, 
as in the normal system. 

In convential superconductors, which contain positively charged ions, paired electrons, and normal electrons at 
T 7^ 0, the photon spectrum is determined by the relation between the vector potential A and the current j, which 
is proportional to A; see Eqs. (96.24) and (97.4) of Ref. |^(|. The analogy with the present case is straightforward. 



"'"Note that in conventional nuclear physics one usually employs particle-hole diagrams even at zero temperature, thereby 
considering nuclear matter to be normal. Small effects of pairing can be neglected, since the typical energy in a nucleonic 
particle-hole diagram is of the order of the Fermi energy ef, and ef > A holds 



3 



From the latter equation, for sufficiently low photon momenta we immediately obtain the relation Airj ~ — m 2 (T)A 
between the Fourier components of the current and the vector potential, where the effective photon mass is 



m 7 (T)^J T<T c , p . (2) 

V m p 

Here m* denotes the effective in-medium proton mass, and p p {T) = p p (T CtP — T)/T c . p denotes the paired proton 
density. The choice of a linear temperature dependence for p p (T) corresponds to the Ginzburg-Landau approach. 

A small complex contribution ~ e 2 f(uj, k)exp(—A p /T)A, where f(oj, k) is a function of the photon frequency u> and 
momentum k, has been neglected in the above relation between j and A. More realistically, for T near T c>p , one 
must take into account this off-shell effect for the photon. At lower temperatures, correction terms are exponentially 
suppressed. Below we take the photon spectrum to be 



lu = y k 2 + m 2 , (3) 

thus neglecting the aforementioned small polarization effects. 

Note that external photons cannot penetrate far into the superconducting region. The photons that we deal with 
are thermal photons with foregoing dispersion law, governed by the corresponding Bose distribution. In consider- 
ing neutrino reactions below, we integrate over the photon phase-space volume thus accurately accounting for the 
distribution of these photons in warm neutron star matter. 

To illustrate more transparently the most important facets of the reconstruction of the photon spectrum in the 
superconducting region, we consider a two-component, locally neutral system consisting of charged fermions (i.e., the 
normal subsystem) described by the Dirac field ip, and a charged condensate (i.e., the superconducting subsystem) 
described by a condensate wave function 

ip = (p c exp (z$). (4) 

The real quantity tp c is the order parameter of the system, i.e., ip 2 ~ n c , where n c is the number density of particles 
in the condensate, and the real value $ is a phase. In a fermionic system with pairing, the density n c is proportional 
to the pairing gap A. 

The equation for the electromagnetic field A p in such a system reads 

□A p = 47r^, (5) 

where the current is 

= eiTpj^ip - ei(ip*d M (p - (fid^ip*) - 2e 2 \(p\ 2 A fl . (6) 

Substituting Eq. (^) into Eq. (j|), the electromagnetic current becomes 

U=tf + 8j», (?) 

where the first term j£ = —2e 2 (p 2 A fJ , is the superconducting current, and the second term Sj p contains the normal 
current j^° T and some response j™ s from the charged condensate, i.e., 

Sj, = + 37 = eih^ + WA$o. (8) 

Due to gauge invariance, the phase 4> = <E>o + < £' is not constrained, and <&o can be chosen in such a way that it cancels 
the normal current, i.e., Sj^ = 0; otherwise the remaining part of the normal current would destroy superconductivity 
and the ground state energy would increase. This compensation of the normal current j" or , which in metals and in 

normal plasma is proportional to the electric field E, is a necessary condition for the existence of superconductivity. 
Only a diamagnetic part of the fermionic current proportional to the electromagnetic field A^ may remain. The 
latter may lead only to a minor (~ e 2 ) contribution to the unit values of dielectric and diamagnetic constants. The 
remaining part of the phase $' is hidden in the gauge field, resulting in the disappearance of the Goldstone field 
(see the analogous discussion of the Higgs effect, e.g., in Ref. |22]). The total number of degrees of freedom does not 
change, so the disappearance of the Goldstone field is compensated by the appearance of an extra (third) polarization 
of the photon. As a result of Eqs. (Q) and (^), the electromagnetic field obeys the equation 

UA p = -8ne 2 ^ 2 c A p , (9) 
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which immediately yields the photon spectrum in the form (tl), where the photon mass is now given by 



m 7 = y/8we 2 (fil. (10) 

What we have demonstrated is known as the Higgs-Meissner effect: in the presence of a superconducting component, 
the photon acquires finite mass. We see that in a two-component (normal + superconducting) system, the photon is 
described by the dispersion relation (||), as it would be in a purely superconducting system, and not by a plasma- like 
dispersion law, as in the absence of superconductivity. Another way to arrive at Eq. (||) is given in the Appendix in a 
non-covariant formulation. Similar derivations for different specific physical systems, guided by the general principle 
of the compensation of the normal currents in a superconductor, can be found in Refs. ||,|l|,|o||| . 

Expressing the amplitude of the condensate field in terms of the paired proton density ||18|, one obtains from 
Eq. ( [lO| ) the result (|J). Taking m*(po) — OSm^ (with mjv the free nucleon mass), with Eqs. (g) and (§) we 

estimate m 7 (/? = po,T) ~ 1.6y/(T cp — T)/T c ^ p MeV <C uj p i(p ~ po). Due to the rather low effective photon mass in 
superconducting neutron star matter at T < T CiP < lo p i , one may expect a corresponding increase in the contribution 
of the above diagrams to neutrino emissivity. 

To avoid misunderstanding, we note the following. At the first glance one might suggest that the photon self-energy 
is completely determined by the above neutrino production diagrams, but with neutrino legs replaced by a photon 
line. If so, the contributions of the electron-loop and proton-loop diagrams would accurately determine the plasmon 
spectrum of photon excitations with energy gap equal to a high plasma frequency (at least if one drops small terms 
proportional to the proton gap in the calculation of the proton-proton-hole diagram, now with an incoming and 
outgoing photon, as suggested for the corresponding neutrino process). How does this relate to the massive photon 
spectrum of superconducting systems? The answer is that in a system with a charged condensate, in addition to the 
cited photon propagation diagrams, there appear specific diagrams for photon rescattering off the condensate given 
by terms oc e 2 </5 2 A p A^, 2e tp\ 9^$ A^ in the corresponding Lagrangian. Their contributions to the equation of motion 
for the electromagnetic field are, respectively, the last two condensate terms in the electromagnetic current in Eq. (|^). 
The specific condensate diagrams responsible for the compensation of the loop diagram contributions in the photon 
propagator make no contribution to neutrino emissivity. Indeed, the neutrino legs cannot be directly connected to 
the photon line via such interactions, (without invoking the internal structure of the condensate order parameter (p c ; 
this contribution is obviously small compared to what we have taken into account). Thus, we have argued that in the 
presence of superconducting protons, neutrino pairs can be produced in the reaction shown by the above diagrams, 
where the photons possess rather small masses generated by the Higgs-Meissner mechanism. 

Having clarified of this important issue, we are ready to calculate the contribution of these processes to neutrino 
emissivity and compare the result with known emission rates. 



III. CALCULATION OF EMISSIVITY 



The matrix element of the above diagrams for the i-th neutrino species (i = {v e , v^, v T }) is 

M (l)a = -iV^e^=e a p (r 7 T«w - l p , (11) 

where 



Tf^-Tr J '*V-T»iG j (p)w} i) '>iG j (p + k), j = {e,p}, (12) 



and 



Gj(p) = (£ + m 3 ) + 2^n» 8{p z - m 2 )6(p ) j (13) 

is the in-medium electron (proton) Green's function, and Uj(jp) = &{pFj — p)- is the corresponding polarization 
four-vector of the massive photon, with three polarization states in superconducting matter. The factor T 7 takes into 
account nucleon-nucleon correlations in the photon vertex. The quantity G = 1.17 -10~ 5 GeV -2 is the Fermi constant 
of the weak interaction. Above, l p denotes the neutrino weak current. The electron and proton weak currents are 

W^o = Y(4 ] - <f 75), = Y{k pp - g Alppl5 ), (14) 
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1.92 and c^ l) = c^ t) = ci ) = 1 - 4 sin 2 6 W 



where c^"' = c ( y ] = 1 + 4 sin 2 6 W 

Integrating Eq. rtlSj) over the energy variable, we obtain for the i-th neutrino species 



0.08. 9w is the Weinberg angle 

= 1. Proton coupling is corrected by nucleon-nucleon correlations, i.e., by the factors k pp and 



• / T w«> _ T Ww>^ = r « + T W pw + r W p/^ 



(15) 



ppp^ ppp 



k 2 [{k-u) 



ppp 



kgu\, 



(16) 



where — (k ■ u)k^ L — u^k 2 , fc M = (ui,k), k 2 — k^k^ = oJ 2 — k 2 . The four-velocity of the medium u M is introduced 
for the sake of covariant notation. The transverse (ft), longitudinal (r;), and axial (t$) components of the tensors in 
Eq. (P) yield 



Ji) 



r£ J = 2c v ) {A e + k 2 B e ) - 2c v >R K (A P + k 2 B p ), 



(-); 



(17) 



r^=r^-r§=Ak^B e -c^R K B p ], 



(18) 
(19) 



where R K — n pp /c v , and 



^"V (2*)»*e0 + 2 1+ 2m 2 r^' 



(20) 



d p rij (p) 



(2tt) 3 



(pfc) 2 



pk 



4E 



(21) 



c? 3 p m 
3 »i(P) 



(27T) 



0)3 



4E, 



W)4 



4*' = / 



(22) 



Here we note that the contribution of the axial component T5 to the resulting neutrino emissivity is small ( T5 jr% ~ 
m 2 T5/(jj 2 Ti ~ m 7 /m^f for protons and ~ (m 1 m e /p 2 Fe ) \n.(pp e /m e ) for electrons), so that it will be omitted. 

The squared matrix element ( |TT| ) for a certain neutrino species, summed over the lepton spins and averaged over 
the three photon polarizations, can be cast in the form 



£l^ (4) l 2 = H G2 



-0)2 



20JlLU2 + 2 



(kqi)(kq 2 
k 2 



(23) 



.0)2 



ujiu>2 + q\<h — 2 



{k-q 1 )(k-q 2 ) „(kqi)(k<&) 



k 2 



- 2- 



k 2 



(24) 



where (k-q\$) = uJU>ifl~ {kqx 2)5 and Wi,2 arid <fi2 denote the frequencies and momenta of the neutrino and antineutrino. 
We have also used the fact that Tr{W"} = 8[q^q'^ + - g^{ qi ■ q 2 ) - Aie^ x Pq lx q2p\- 
The emissivity of our processes is given by 



G 



d 3 k 



9i 



rf 3 



<h 



(2tt) 3 2lu (2tt) 3 2wi (27r) 3 2w 2 / ^1+^2 



(25) 



r 



- 1 



£ ^|MW| 2 (27r) 4 5 4 (fc- gi -g 2 ). 



(26) 



Substituting Eq. (||) into Eq. (Kg), we finally obtain 



e2 = _,T° 7re 2 G 2 a 2 J, / = 



9(2tt) 



ef-1 



..(f?) + 'f(fa) 



(27) 



7TL 

where a = -j?-, and 



y K 2m; 1 + ] v [ 8ir Pp) [ 2> 



(28) 



r 2 (,)«4x 2 J] 



2ml 



(29) 



Some numerically small terms have been dropped in Eq. (j28|). 

The integral I can be calculated analytically in the two limiting cases a -C 1 and a S> 1: 



/(a > 1) 



|E Q 3/a ( i + |. )e -« [r 2 (1)+r 2 (1)]) 



(30) 



/(a « 1) « 2C(3) h(0) + r 2 (0)], C(3) ^ 1.202. (31) 



Thus, combining Eqs. (|l|) and (|27|)-(|30|), we obtain an estimate for the emissivity of our reactions (we present here 
the result for m 7 > T and for three neutrino species): 



el « 2.6 • 10- Tl e-^(^Y{ JL ) 7 1 + [1 + ,] (32) 

VJVleV/ \Po/ \ 2m 7 / cm J sec 
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„ = 0.0003^ ( = ) [j-j -0.035*^^-] . ,33, 



Here Tg denotes temperature measured in units of 10 9 K. The 1 in square brackets in Eq. (32) corresponds to the 
electron-electron- hole diagram, whereas the factor 7/ is related to the proton-proton-hole (first term in Eq. j33|)) and 
the interference diagrams (second term in Eq. (j33|)). 



The emissivity Eq. (32) varies with temperature as T 3 / 2 exp(— m 7 /T), whereas the emissivity of the modified URCA 
process varies as T s exp(— (A p +A„)/T) in the region of proton (A p ^ 0) and neutron (A n ^ 0) pairing. Hence, one can 
expect that the process j m — * vv will dominate at comparatively low temperatures, when A p (T) + A n (T) — m~ ( (T) > 
and T < T CtP . 

IV. NUMERICAL ESTIMATES 

To obtain quantitative estimates we need the values of the nucleon-nucleon correlation factors k pp and T 7 . According 
to Ref. |Q we can exploit 

k pp = Cy ^ 2 f n pG\)A nn T (f nn ) , (34) 
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where f np ~ —0.75 and f nn ~ 1.25 are the constants in the theory of finite Fermi systems (T^^^]; C 1 = m^pFn/^ 2 
is the density of states at the Fermi surface; A nn is the neutron-neutron-hole loop, 

Co A nn = zCo / G„(p + k) G n ( P ) « (35) 



for values of w » |fc|j)ir„/m* of interest, and T 1 (f n n) = 1 — 2f nn CoA 



nn • 



We note that the second term in Eq. ([34|) is not proportional to a small factor c v \ because the nucleon-nucleon 
correlations also allow for the emission of j/f'-pairs from the nnT x loop. Numerical estimates of the ratio R K are 
as follows: for a ^> 1, R K ~ 1.6 for p = pq, m*(po) — 0.8m n , and R K ~ 2.1 for p = 2/?o, m*(2po) — 0.7m„; for 
a <C 1, -R K — 1 and correlation effects are negligible. The in-medium renormalization of the proton electric charge 
included in the factor T 7 can be also expressed in terms of the constants in the theory of finite Fermi systems and 
the proton-proton loop factor (A pp ); see Ref. jF9|. The latter is suppressed at relatively low proton densities. We 
can therefore take T 7 w 1. With these estimates, we observe that the main contribution to neutrino emissivity comes 
from electron-electron-hole processes. 

The ratio of the emissivity el (|32| ) to the emissivity e^ M of the modified URCA process, Rfm = el/ £ v M ■ 1S 

3 

A 



Rfm « 1.5- 10 4 -T ( 



For further estimates we need the values of the neutron and proton gaps, which are unfortunately model-dependent. 
For instance, the evaluation in Ref. Q yields A„(0) ~ 8.4T Ci „ ~ 0.6 MeV, T c . n ~ 0.07 MeV for 3P 2 neutron 
pairing at p = p , and A p (0) ~ 1.76 T CiP ~ 3 MeV, T c ^ p ~ 1.7 MeV for 15 proton pairing, while Ref. |27| uses 
A n (0) ~ 2.1 MeV, T Cj „ ~ 0.25 MeV and A p (0) ~ 0.7 MeV, T CjP ~ 0.4 MeV for p = p . Employing these estimates 
of the zero-temperature gaps, their temperature dependence, and the photon effective mass, we obtain from Eq. ( |3^ ) 
the temperature dependence of the ratio Rfm- 

In order to find the lower temperature limit at which the processes j m — > vv are still operative, we need to compare 
the value e~l with the photon emissivity at the neutron star surface, e 7 = 3aT^/R, where a is the Stefan— Boltzmann 
constant, T s denotes the surface temperature of the star, and R is the star's radius. By employing a relation pS[ ] 
between the surface and interior temperatures, we obtain for R 1 = e1js a 



R 1 w 1.2 • 10 9 r 9 -°- 7 e T 



where the star's radius and mass are taken to be 10 km and 1.4M©, with M the solar mass and p some averaged 
value of the density in the neutron star interior. 

The ratios Rfm and i? 7 are plotted as a function of the temperature in Fig. 1 for both of the foregoing parameter 
choices. We see that our new processes are operative in the temperature range 1 • 10 9 K< T < 8 • 10 9 K for the 
parameter choice of Ref. p6]|, and 1 • 10 9 K< T < 4 • 10 9 K for the parameters of Ref. p^J. As one observes in Fig. 1, 
within these intervals the new cooling channel might exceed known cooling processes by up to a factor 10 6 . 

V. CONCLUDING REMARKS 

As mentioned above, for T > T Cj „,T CjP , i.e., in a normal plasma region of the star's crust and interior, photons 
with approximately the electron plasma frequency^] u> p \ can decay into neutrino pairs, as has been shown in previous 
estimates JToj ] . At T < T ClP , however, we are already dealing with massive photons in the region of proton pairing, 
and our new reaction channels can significantly contribute to cooling. 

Our processes can also occur in a charged-pion (or kaon) condensate state but they are suppressed due to the high 
effective photon massp*| m 7 ~ ^8ne 2 (pl ~ 6 MeV for the condensate field <p c ~ 0.1m ff ~ 14 MeV. 

In deriving the value of e^ M used above, one describes the nucleon-nucleon interaction essentially by free one-pion 
exchange. In reality, however, at p > (0.5 — l)po the total nucleon-nucleon interaction does not reduce to free one- 
pion exchange, because of the strong polarization of the medium, whereby a significant part comes from in-medium 



§ A rather small extra contribution also comes from the proton-proton-hole diagram. 

"For simplicity, in this estimate the peculiarities of a condensate with nonvanishing momentum M are ignored. 
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pionic excitations J7|,p4|]5||| . Occurring in intermediate states of the reaction, the in-medium pions can also decay 
into ev, or first into a nucleon-nucleon-hole, which then radiates ev, thereby substantially increasing the resulting 
emissivity. Other reaction channels such as n — > n pa i r vv and p — > p pa irvv open up in the superfluid phase with paired 
nucleons [^5| , [?4| ](t| . All these reaction channels give rise to a larger contribution to the emissivity than that of the 
modified URCA process estimated via free one-pion exchange. Above we compared el with e^ M just because the 
latter is used in the standard scenarios of neutron star cooling. 

As we also mentioned in the Introduction, there are other processes like those considered above. Emissivity of the 
process pj m —> P P airVv is substantially suppressed (at least by a factor e 2 and also due to a much smaller phase-space 
volume) compared to that of the process p — > p pa irvv. According to simple estimates, e.g., using Eq. (22) of Ref. [[l6| , 
the process cy — > evv makes a very small contribution to the emissivity both in the inner crust and in the interior 
of neutron stars, even when one neglects the photon mass. Thus we may conclude that the process ej m — » evv also 
leads to a minor contribution to the emissivity at the densities and temperatures under consideration. 

In summary, the processes j m — + ee _1 + pp^ 1 — ► vv might be operative over some temperature interval 
T ~ (10 9 — 10 10 ) K, T < T CiP , and together with other in-medium modified processes ||, they should be incor- 
porated into computer simulations of neutron star cooling. 
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Science Foundation and the Russian Government. B. K. and E. E. K. are supported by BMBF Grant 06DR666. 
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APPENDIX 



We can also achieve the same results that led to Eq. (10) by starting with Maxwell's equations (in obvious notation): 



ik- E = Airp, ik x B = 4vrJ - itoE, (38) 



k-B = 0, kxE = uB, (39) 

where the charge density p is the superposition of the density of free charges and the density of bound charge. Full free 
charge density being zero in our case due to local electro-neutrality. The current j is a superposition of an external 
test current and the induced current 

J = J cxt +J md - 

In normal systems, the induced current (i.e., the current of non-paired charged particles) j md = j noT is related 
to E via longitudinal and transverse dielectric constants e/ and e t . This connection results in longitudinal and 
transverse branches of the electromagnetic excitations, with an effective photon gap equal to the plasma frequency 
uj p i [ jTo| . In contrast, in a superconducting system the condensate makes two other contributions to the current, namely 

j A = —2e 2 (p 2 A and j res = 2e<p 2 V<I>. Letting $ = $(1) + < &(2) , we have j rcs = im +jT|)- These two terms are determined 
as follows. As we have argued above, superconductivity requires the compensation of the normal component of the 
current proportional to E, i.e., we can take j nor + jffi = 0. Only small contributions ~ e 2 exp(— A p /T)lu 2 A and 

~ e 2 exp(— A p /T)k 2 A, as well as a small imaginary contribution ~ ie 2 F(ui, k)ex-p(—A p /T)A, where F is some function 
of oj and k, can still remain from the value j nor (see Eqs (96.24) and (97.4) of Ref. |p0|] ). We neglect these small 
contributions. The part of the current ~ V<5(2) can be hidden in j A by a gauge transformation of the field A. We 
then have 

ik x B ~ j A - iojE. 

Taking the vector product of this equation with k, we obtain (ui 2 — k 2 — 8ire 2 tp 2 )B = 0. From this relation we observe 
that the electromagnetic excitations possess the mass given by Eq. @. Hence, we have demonstrated that one can 
obtain the well-known plasma photon spectrum for a normal system, and at the same time one can obtain a massive 
photon spectrum and the Higgs-Meissner effect in a system with a charged condensate. 
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FIG. 1. Temperature dependence of the ratios Rfm and _R 7 at nucleon density p = po- Solid curves correspond to the 
parameter choice of Ref. [^6| , whereas the dashed curves depict results with parameters of ref. . Shaded bars indicate the 
temperature regions in which cooling via massive photon decay is more efficient than standard cooling processes. 
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